Fitting to data of superluminal neutrinos with phenomenological scenarios 
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We test several phenomenological secnarios of faster-than-light neutrinos, by fitting to 
the experimental data from OPERA, MINOS and Fermilab79. Om- purpose is to see, from 
the perspective of the current data, whether or not the speed of the superluminal neutrino 
depends on its energy. We show that the Coleman- Glashow scenario in which the velocity of 
the neutrino is free of the energy fits the data best. However, the result of SN1987A cannot 
be explained by this model. We find that a power-law model with the power close to zero 
can simultaneously explain the results of SN1987A and OPERA+MINOS+Fermilab79. 

The OPERA collaboration reported the superluminal muon neutrino {u^) data [1], recently. 
The neutrinos with average energy ~ 17 GeV arrive earlier by 

6t = (60.7 ± 6.9 (stat.) ± 7.4 (sys.)) ns (1) 

than photons, from CERN to Gran Sasso Laboratory with distance about 730 km. This indicates 
that neutrinos are superluminal, with an excess of the speed than light (in vacuum) 

v-l = (2.48 ± 0.28 (stat.) ± 0.30 (sys.)) x 10"^ (2) 

with significance level of Qa. Note that throughout the paper we use the natural units with 
c = 1. It was also found by OPERA that the velocity difference — 1 is almost independent of 
energy, by dividing the events into two groups with energies below or above 20 GeV: v — 1 = 
(2.16 ± 0.76 (stat.) ± 0.30 (sys.)) for {E) = 13.9 GeV and i; - 1 = (2.74 ± 0.74 (stat.) ± 0.30 (sys.)) 
for {E) = 42.9 GeV. 

It should also be pointed out that earlier experiments ever obtained similar results of superlu- 
minal neutrinos, though with lower significance. For example, the MINOS experiment at Fermilab 
in 2007 [2] and even earlier experiments at Fermilab in 1979 [3] reported the results of f — 1 ~ 10~^, 
similar to that of OPERA. These results are summarized in Table I. However, the supernova neu- 
trinos from SN1987A place a stringent velocity constraint, u — 1 < 2 x 10~^, for electron neutrinos 
with energies from 5 to 40 MeV, seemingly inconsistent with the superluminal neutrino result of 
OPERA. 

At first glance, the result of that the neutrino travels faster than light is too shocking, and 
many people instinctively reject such a result. Indeed, there might be some systematic errors 
or some other unknown factors in the experiment that have not been taken into account. For 
example, it was suspected that there may be some serious problems in the measurement of time 
and distance in the experiment [4]. And, in the theoretical respect, Cohen and Glashow [5] argued 
that superluminal neutrinos with high energies would lose energy rapidly via bremsstrahlung of 
electron-positron pairs, causing the beam to be depleted of higher energy neutrinos, and so they 
refuted the superluminal interpretation of the OPERA result. However, on the other hand, if the 
nature of neutrinos is indeed superluminal, then the modern physics would be revised enormously. 
Thus, before the independent experiments, such as BOREXINO, ICARUS, and MINOS, confirm 
(or negate) the OPERA result, it is a rational way for us to suppose that the data and conclusion 
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TABLE I: The summary of supcrluminal neutrinos from OPERA, MINOS, and Fermilab79. 



Experiment 


Velocity constraint 


Energy range 


Reference 


OPERA 


(2.48 ± 0.58) X 10"^ (6.0 cr) 


- 17 GeV 


[1] 


MINOS 


w - 1 (5.1 ± 2.9) X 10-^ (68% CL) 


- 3 GeV 


[2] 


Fcrmilab79 


|w - 1| < 4 X 10-^ (95% CL) 


30 to 200 GeV 


[3] 



of OPERA are right and try to look for the possible clues of the new revolution of physics via the 
super luminal neutrino experiments. 

In recent days, many theoretical explanations for the superluminal neutrinos have been put 
forward. For example, it has been suggested that such a phenomenon might originate from a 
violation of Lorentz invariance [6-9]. Nevertheless, it is hard to understand why the Lorentz 
violation happens at such low energies, far below the energy scale of quantum gravity. Tachyonic 
neutrino scenario is another attractive explanation. However, this scenario is clearly disfavored by 
the data [6]. Besides, there are also many other possible mechanisms, such as, extra-dimension 
model [10], neutrino dark energy model [11], and so on. 

One important point should be clarified, in our opinion, before we make deeper theoretical 
analysis: whether or not the superluminal phenomenon of neutrinos is energy-dependent? In fact, 
it was proposed by Li and Wang [7] that the Lorentz violation might be energy-dependent, only 
occurring in a window of particular energies (including 3 — 200 GeV) , thus reconciling the results of 
OPERA and SN1987A. In addition, Amelino-Camelia et al. [6] also tackled the issue of the energy- 
dependence of the speed of neutrinos. In [6] several cases of violation or deformation of Lorentz 
invariance were considered, and the authors showed that the Coleman-Glashow (CG) scenario and 
the Doubly Special Relativity (DSR) scenario with linear-dependence are more favored, while the 
tachyon scenario and DSR model with quadratic-dependence are clearly disfavored, by the data. 
Of course, it seems that the data-favored models are inconsistent with the result of SN1987A. In 
this paper, we will make a more detailed analysis of current data with some phenomenological 
scenarios. We will show whether the superluminal phenomenon of neutrinos is energy-dependent, 
according to the current experimental data, and which phenomenological scenario can fit the data 
well and be consistent with the result of SN1987A at the same time. 

We show the data of OPERA+MINOS+Fermilab79 in the left panel of Fig. 1. As discussed 
in [3, 6], when using the data of Fermilab79, a bias correction, downward shift of data with 
61979 = 0.5 X 10~^, should be considered. The data with the above correction (for Fermilab79) are 
shown in the right panel of Fig. 1. 

In the ultra-relativistic regime, the Einstein's mass-energy relation gives the following descrip- 
tion of the dependence of speed on energy: 

1 , , 

Obviously, in the theory of relativity a normal particle with real mass m can never travel faster than 
light in the vacuum. Now, let us estimate the order of magnitude of the velocity difference for the 
neutrino by using Eq. (3). If the mass of neutrino is supposed to be m ~ eV and the neutrino energy 
E under consideration ranges from 3 GeV to 200 GeV, we obtain m?/2E'^ ~ 10~^* — 10~^^. So, 
even if we consider a tachyon model in which the neutrino has an imaginary mass, = —rn?, the 
theoretical value obtained is many order of magnitude less than the experimental result, O(10~^). 

Supposing that there is a theoretical mechanism that can explain the experiments such as 
OPERA, the above relation can be modified, phenomenologically, as 

v{E)-l = -\'^ + f{E)^f{E), (4) 
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FIG. 1: The superluminal neutrino data from OPERA, MINOS and Fermilab79. The left panel shows the 
raw data, and the right panel shows the data corrected by 61979 = 0.5 x 10"'^ for Fermilab79. 



where the form of f{E) is determined by the underlying theoretical mechanism that is unknown 
for us now but can be modeled phenomenologically. We can view that f{E) comes from a Lorentz 
violation, though we cannot understand why the energy scale of such a Lorentz violation is so 
low. Alternatively, we can even consider a mass-running tachyonic neutrino scenario with mass 
^ _m{Ef = 2E'^f{E). 

In what follows, we take a totally phenomenological perspective, assuming the possible forms 
of f{E). First, we consider the possible form inspired by Coleman and Glashow [12, 13], 

f{E) = 6, (5) 

where 5 is a constant parameter. Also, we follow Amelino-Camelia et al. [6], taking the DSR- 
type cases, f{E) = iiE and f{E) = i2E'^. We, furthermore, consider a more general form, a 
combination of the above possible cases, 

f{E) = aE'^ +^E + S. (6) 

We will use the superluminal neutrino data currently available, namely, the combination of OPERA, 
MINOS and Fermilab79, to tell us which scenario works well. To fit the data, we use a statistic, 

^^^^ = 1^ ' 

1=1 

where cr{Ei) is the la error for each datum of v{Ei), and denotes the model parameters. By min- 
imizing x^, we can find the best-fitted parameters of the models, and further obtain the probability 
contours in the parameter-planes. 

For the CG scenario, we get xLii = 56.94. The DSR-type cases, f{E) = hE and f{E) = i^E"^, 
are both clearly disfavored by the data, leading to much bigger Xmin' being 195.46 and 368.18, 
respectively. The combination form, f{E) = aE"^ + f3E + 5, gives the least Xmin' 56.77, among these 
scenarios. However, it is unwise to use the statistic alone to compare competing models since 
the number of parameters is different for the models. In general, a model with more parameters 
tends to give a lower Xmin' instead one may employ the information criteria (IC) to assess 
different models. In this paper, we use the BIG (Bayesian information criterion) and AIG (Akaike 
information criterion) as model selection criteria, defined respectively as 

BIC = Xmin + A:ln7V, AIC = xL„ + 2^, (8) 
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TABLE II: The fit results of the phcnomcnological scenarios. 



f(E) scenario 


Amm 


k 


ABIC 


AAIC 


5 


56.94 


1 








liE 


195.46 


1 


138.52 


138.52 




368.18 


1 


311.24 


311.24 


aE'^ +PE + 6 


56.77 


3 


4.62 


3.83 


^{E/GcVy 


56.78 


2 


2.24 


1.84 



where k is the number of parameters, and is the data points used in the fit. Note that the 
absolute value of the criterion is not of interest, only the relative value between different models, 
ABIC or AAIC, is useful. For the details about the BIC and AIC, especially their applications 
in a cosmological context, see, e.g., [14, 15]. According to the fitting result, we see that the CG 
scenario has the minimal values of BIC and AIC. The scenario of f{E) = aE"^ + I5E + 6 yields 
ABIC = 4.62 and AAIC = 3.83, relative to the CG case; see Table II. So, according to the 
principle of Occams razor, "entities must not be multiplied beyond necessity," the CG scenario, 
f{E) = 5, is preferred by the data. Figure 2 shows the fit result of the CG scenario, with the 
left panel the one-dimensional likelihood distribution of the parameter 6, and the right panel the 
best-fitted case of v{E) comparing to the data. Figure 3 shows the fit result of the scenario of 
f(E) = aE'^ + f3E + 5. From this figure, we can see that, though the parameters a and /3 are both 
around zero at the best fit, there are some evident degeneracies between the parameters. From the 
probability contours in the parameter-planes, we find that, a and f3 are anti-correlated, /3 and 6 
are also anti-correlated, and so a and 6 are in positive correlation. The best-fitted v{E) curve is 
similar to a horizontal line, but with a slight slope in the right-hand region. 




FIG. 2: Likelihood distribution of parameter S for the Coleman-Glashow scenario, f{E) = S. 

While the CG scenario can fit the superluminal neutrino data well, it cannot simultaneously 
explain the result of SN1987A. Can we find a phenomenological scenario that can reconcile the 
SN1987A with the superluminal neutrino experiments such as OPERA, MINOS and Fermilab79? 
For this purpose, we test the following power-law scenario, 

fiE) = CE\ (9) 

where the energy E is in units of GeV. We now fit to the experimental data of OPERA, MINOS 
and Fermilab79 with this model. From Table II we see that the power-law model yields a lower 



5 




-0003 0.000 0.003 20 40 60 80 100 120 140 160 180 200 

p (Gev"') E(GeV) 

FIG. 3: Probability contours at 68.3%, 95.4% and 99.7% confidence level in a — /3, a — (5 and (5 — 5 planes 
for the scenario, f{E) — aE^ + j5E + 5. 

^min' 56.78, than that of the CG model. However, when considering the number of parameters, 
the power-law model produces slightly larger BIG and AIG values than the GG model, namely, 
ABIG = 2.24 and AAIG = 1.84. This indicates that the power-law scenario can also fit the data 
very well. We show the probability contours in ^ — A plane in the left panel of Fig. 4. It can be seen 
from this figure that the best-fitted A is a tiny positive number, in the vicinity of zero. In addition, 
the parameters ^ and A are in evident degeneracy, strongly anti-correlated. The right panel of 
Fig. 4 shows the best-fitted v{E) curve of the power-law scenario. We can see that this scenario 
behaves just like the GG scenario, and the only difference lies in the region of energy approaching 
zero, i.e., in the power-law model the speed of neutrino will rapidly decrease to the speed of light 
as the energy goes to zero, and so the result of SN1987A can also be successfully accommodated 
in this scenario. 

In summary, we have tested several phenomenological scenarios of superluminal neutrinos, ac- 
cording to the fits of the experimental data of OPERA, MINOS and Fermilab79. The purpose 
of this work is to see whether or not the speed of the faster-than-light neutrino is dependent on 
the energy, from an analysis on the current data. We adopted a phenomenological perspective 
that the only focus is on the data but not the deeper theoretical mechanism. We have shown 
that the GG scenario in which the velocity of neutrino is free of the energy can fit the data best. 
Nevertheless, this model cannot meanwhile provide an explanation for the result of SN1987A. We 
found that a power-law scenario with the power being a tiny positive number can simultaneously 
explain the results of SN1987A and OPERA -|-MINOS-|-Fermilab79. The scenario may deserve a 
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FIG. 4: Probability contours at 68.3%, 95.4% and 99.7% confidence level in ^ — A plane for the power-law 
scenario, f{E) = ^E^. 



further investigation. 



Acknowledgments 

We thank Nan Li, Yang- Yang Li, Yun-He Li and Jing-Fei Zhang for useful discussions. This work 
was supported by the National Science Foundation of China under Grant Nos. 10705041, 10975032, 
and 11175042, and by the National Ministry of Education of China under Grant Nos. NCET-09- 
0276 and N100505001. 




[1] T. Adam et al. [OPERA Collaboration], arXiv:1109.4897 [hep-cx]. 

[2] P. Adamson et al. [MINOS Collaboration], Phys. Rev. D 76, 072005 (2007) [arXiv:0706.0437 [hep-ex]]. 

[3] G. R. Kalbfleisch et al. Phys. Rev. Lett. 43, 1361 (1979). 

[4] R. A. J. van Elburg, arXiv:1110.2685 [physics. gen-ph]. 

[5] A. G. Cohen and S. L. Glashow, arXiv:1109.6562 [hep-ph]. 

[6] G. Anielino-Camelia, G. Gubitosi, N. Loret, F. Mercati, G. Rosati and P. Lipari, arXiv:1109.5172 
[hep-ph]. 

[7] M. Li and T. Wang, arXiv: 1109.5924 [hcp-ph]. 

[8] Z. LingU and B. Q. Ma, arXiv:1109.6097 [hep-ph]. 

[9] X. J. Bi, P. F. Yin, Z. H. Yu and Q. Yuan, arXiv:1109.6667 [hep-ph]. 

[10] S. S. Gubscr, arXiv:1109.5687 [hep-th]. 

[11] E. Ciuffoh, J. Evslin, J. Liu and X. M. Zhang, arXiv:1109.6641 [hep-ph]. 

[12] S. R. Coleman and S. L. Glashow, Phys. Lett. B 405, 249 (1997) [arXiv:hep-ph/9703240]. 

[13] S. R. Coleman and S. L. Glashow, Phys. Rev. D 59, 116008 (1999) [arXiv:hep-ph/9812418]. 

[14] A. R. Liddle, Mon. Not. Roy Astron. Soc. 351, L49 (2004) [arXiv:astro-ph/0401198]. 

[15] M. Li, X. D. Li and X. Zhang, Sci. China Phys. Mech. Astron. 53, 1631 (2010) [arXiv:0912.3988 
[astro-ph.CO]]. 



